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We consider a multipartite system consisting of two noninteracting qubits each embedded 
in a single-mode leaky cavity, in turn connected to an external bosonic reservoir. Initially, 
we take the two qubits in an entangled state while the cavities and the reservoirs have 
zero photons. We investigate, in this six-partite quantum system, the transfer of quantum 
discord from the qubits to the cavities and reservoirs. We show that this transfer occurs 
also when the cavities are not entangled. Moreover, we discuss how quantum discord can 
be extracted from the cavities and transferred to distant systems by traveling leaking 
photons, using the input-output theory. 

Keywords: Multipartite open quantum systems; dynamics of quantum correlations; ex- 
traction of quantum correlations. 



1. Introduction 

In the study of correlations present in quantum states two different approaches 
are nowadays typically distinguished: one is based on the entanglement-versus- 
separability paradigm and the other on a quantum-classical dichotomy.^ In the 
entanglement-separability paradigm introduced by WerneJ^I, the state of a multi- 
partite system is named separable if it can be represented as a mixture convex 
combination of product states relative to the various parts of the total system, oth- 
erwise, it is termed entangled. Entanglement is considered a key ingredient in the 
increased efficiency of quantum computing compared to classical computation for 
certain quantum alg orit hms and thus it plays a central role in quantum information 
and communication.lSEl 

However, entanglement does not exhaust the realm of quantum correlations. A 
quantum state of a composed system may contain other types of nonclassical corre- 
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lation even if it is separableP In the framework of quantum-classical dichotomy, the 
total correlations present in the system can be separated in a purely quantum part 
and a classical part. A measure of quantum correlations is the quantum discord I^'^ 
A mixed bipartite separable state with no nzero quantum discord may yet be ex- 
ploited in quantum computation protocolsl^^SEl intense research activity has 
been dedicat ed to the characterizatio n of quantum discord for several classes of 
both bipartitPHmH and rnu ltipartit quantum states, even in the case of con- 
tinuous variable 

systemJEHSl. 

An experimental observation of the quantu m d iscord 
in nuclear magnetic resonance quadrupolar system was recently reportedU^ 

The dynamics of entanglement for bipartite quantum systems interacting with 
independent or common environments, cither Markov ian or non-M arkovian, with 
the appearance of phenomena like sudden deathP^^U^, revivalJI^I^ or trappinj^, 
can be considered well understood in its general lines. Dynamics of quantum dis- 
cord has also r eceive d much attention f or two-q ubit systems in the presence of 
both Markoviad22112| 

and non- 

Markoviail23124i25| 

environments. Interestingly, dif- 
ferently from entanglement, Markovian evolution can never lead to a sudden death 
of discord. Quantum discord docs result in fact to be more robust than the entan- 
glement against dccoherence. In some cases, in front of a decay of total correlations, 
classical and quantum correlations may present finite time intervals when each of 
them remains constant .l^SI 

Recently, the transfer of entanglement b etween th e parties of a composed system 
has been analyzed in several investigations.^^ZEHEH When two initially correlated 
atoms are placed in two noninteracting leaky cavities, each connected to its own 
reservoir, it has been fo und that entanglement can be transferred from atoms to 
reservoirs via the cavities.^^In this case, different regimes have been shown to exist 
and in particular, under certain conditions, the cavities do not become entangled 
during the dynamics. The reverse problem of entanglement transfer from radiation 
modes to qubits through cavities has been also investigated, findin g the conditions 
for exchange of quantum correlations among the subsystems .^^^^^ Here, we investi- 
gate the conditions for transfer of quantum correlations from qubits inside cavities 
to outside systems. In view of exploiting quantum discord, in this paper we analyze 
how it can be extracted for this kind of systems. The characteristics of discord 
transfer shall be compared with the ones known for entanglement transfer. 



2. Model 

We consider a system composed by two noninteracting subsystems {i ~ 1,2), each 
consisting in a qubit (two-level emitter) qi coupled to a single-mode cavity Ci in 
turn interacting with an external reservoir Vi (see Fig. [1]) . The Hamiltonian of the 
total system is thus given by the sum of the Hamiltonians of the two noninteracting 
subsystems 
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subsystem 1 subsystem 2 



Fig. 1. (Color online) Schematic representation of the six-partite system. The two qubits qi and 
92 arc initially entangled. 

In each subsystem, we distinguish the two-partite system of interest S made by the 
qubit plus the cavity from the reservoir r made by the external bosonic modes. The 
Hamiltonian of each part i = 1, 2 reads like (we omit index i, h ~ 1) 

H = Hs + Hr + Hi , Hs ^ ^^^of^z + a^a + g{a^ + (7+ a) , 
Hr = Y,^khlbk, Hi ^Y.^'kiabl + aHk), (2) 

k k 

where g is the coupling constant between qubit and cavity, the coupling constants 
between cavity and external modes, CTz denotes the usual diagonal Pauli matrix, 
cr± arc the two-level raising and lowering operators, a and bk are the annihilation 
operators for the cavity and for the reservoir modes. 

2.1. Dynamics of subsystems 

Being the two subsystems noninteracting, they evolve independently so that we 
can analyze the dynamics of only one subsystem and use it to obtain the evolution 
of the global six-partite system. We will consider initial states such that, in each 
subsystem, only a single excitation is initially present in the qubit. This means that 
we will need to know the evolution of the single subsystem state \(po) — |l)(j|0)c|O)r, 
where |0)r = 11^10^)^. Under the action of the Hamiltonian of Eq. ([2]), \ipo) evolves 
as 

IVt) = (t\l)q\0)c\O)r + 77t|0>,|l)c|O), + SfeAfc(t)|G),|G),|lfe)r, (3) 

which can be rewritten in terms of a collective state of the reservoir modes. \l)r — 
{l/Xt)^kMt)\U)r, as 

l^t) = 6|l)g|0)c|O). + r;t|0),|l),|O), -f- xt|0),|0),|l),. (4) 

Eq. allows to compute the joint evolution of the total six-partite system starting 
from an arbitrary initial state where only one excitation is present in each atom, 
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as we will show explicitly in Sec. El In the following we obtain the time dependent 
coefficients ^t, rjt and Xt, which may in general be different for the two subsystems. 
For each subsystem |xtP = 1 — |^tP — |?7tP, so that it is sufficient to have and 
rjt- To this aim, we need to solve the reduced dynamics of the bipartite system of 
interest S (qubit plus cavity) and compare its density matrix at time t with the 
one obtained by tracing out the reservoir in the evolved state of Eq. (j4|), whose 
elements in the qubit-cavity basis B {|1) = |11), |2) = |10), |3) = |01), |4) = |00)} 
are 

PS22 = PS33 = I'ytP, PS44 = IXtP, PS23 = PS*32 = 6'7t ■ (5) 

The reduced density matrix at time t of the qubit-cavity system S, induced by the 
Hamiltonian of Eq. ([2|), may be obtained using a phenomenological master equation 
of the form 

^PS = i[ps,Hs] + '^{2apsa'^ - a)aps ~ Pso^a), (6) 

where 7 represents the rate of loss of photons from the cavity. The above master 
equation is appropriate when the reservoir is at zero temperature, the coupling 
between the cavity and the external modes of the reservoir has a flat speetrumin 
the range of involved frequencies and the qubit is resonant with the cavity I 33 | 34 | 
We will limit our investigation to this physical condition, adopting therefore the 
master equation of Eq. ([6]). The comparison between the solution of this master 
equation with the initial condition |l)g|0)c and Eq. ([5]) allows to obtain and T}t 
as 



1^,1^ =e-* cos(m) + ^ sin(r!t) , |r,,f^ = ^e"* sin^(m), (7) 



AO. 



^ 9' 



where we have introduced the characteristic frequency = \/ ~ 7-^/16. When 
(4g)/7 > 1 the functions |^t| and |77t| have a damped oscillatory behavior, while 
when (4(7)/7 < 1 they become hyperbolic and oscillations disappear. When qubit 
and cavity are slightly out of resonance, the master equation of Eq. (|6l) is expected 
to hold as well but the expressions of and r\t become quite cumbersome. 

In order to investigate the flow of classical and quantum correlations within the 
multipartite system, we first introduce some quantifiers able to distinguish quantum 
from classical correlations. 



3. Classical and quantum correlations 

As said before, entanglement does not necessarily exhaust all quantum correlations 
present in a state. An attempt to quantify all the nonclassical correlations present 
in a system, besides entanglement, has led to the introduction of the quantum 
discord, given by the difference between two expressions of mutual information 
extended from classical to quantum system^^^. Following this framework, the total 
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correlations between two quantum systems A and B are quantified by the quantum 
mutual information 

lipAB) ^ S{pa) + S{pb) - S{pab) (8) 

where S{p) = — Tr(plog2 p) is the Von Neumann entropy and pa(b) = ^'^b(A){pab)- 
It is largely accepted that quantum mutual information I{pAB ) is the information- 
theoretic measure of total correlations in a bipartite quantum stati^B] 

On the other hand, the classical part of correlations is defined as the maximum 
information about one subsystem that can be obtained by performing a measure- 
ment on the other system. Given a set of projective (von Neumann) measurements 
described by a complete set of orthogonal projectors {11^} and locally performed 
only on system B, the information about A is the difference between the initial en- 
tropy of A and the conditional entropy, that is2{pAB\{^k}) = S{pA) — ^kPkS{pk), 
where pk = {1 ® 11^)^(1 ® nfc)/Tr[(I ® 11^)^(1 ® 11^)], pk is the probability of 
the measurement outcome k and I is the identity operator for subsystem A. 
To ensure that one captures all classical correlations, one needs to maximize 
1-{pAB\{^k]) over all the sets {11^}. Classical correlations are thus quantified by 
Q{pab) = sup|]-[^jT(pyiB|{nfc}) and the quantum discord is then defined by 

V{pab) = Apab) - Q{pab), (9) 

which is zero only for states with classical correlations and nonzero for states with 
quantum correlations. The nonclassical correlations captured by the quantum dis- 
cord may be present even in separable states.^ 

The maximization procedure involved in computing the quantum discord has 
been analytically solved for certain class of quantum states. Along the paper, we 
will limit our treatment to a subclass of X-structured density operators which, 
in the standard two-qubit computational basis B = {|1) = |11),|2) = |10),|3) = 
|01), |4) = |00)}, has density matrix elements given by 

Pii = a, P22 = /533 = b, Pa = d, pu = w, p23 = z, (10) 

where the coherences pi4 and P23 are real numbers and P22 = P33, so that 
S{pa) = S{pb)- Under this condition, the classical correlations assume the same 
value irrespective of wh ethe r the maximization procedure involves measureme nts 
performed on A or on B.^Elpor this class of states, quantum discord is given bj^^ll 

V{p) = mui{Di,D2}, 



a 



Di = S{pa) - S{pab) - alog2 {^-^—f^ ] - blog^ 
~^1°S2 (7-7-7 )- ^log2 ^ 



D2 = S{pa) - SipAB) - A+ log2 A+ - A+ - log2 A_, (11) 

with A± = i(l ± r) and F^ = (a - d)^ -1- 4(|z| + \w\)'^. Using this expression for 
quantum discord, classical correlations can be in turn obtained by Eqs. © and (|5]). 
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In our anal ysis wc will compare quantum discord and entanglement, quantified 

Toy] 

by concurrence^ which varies from C = for a disentangled state to C = 1 for a 
maximally entangled state. For states defined as in Eq. (fTO|) . concurrence results to 
be C = 2max{0, z - V^, w - b}^ 

4. Correlations transfer by input-output theory 

The dynamics of entanglement between two noninteracting qubits, each in its own 
cavity, has been analyzed by exploiting a model where each qubit is cou pled to a 
continuous ensemble of cavity modes with a Lorentzian spectrumff^I^Hl is pos- 
sible to show that this model is formally equivalent to the quasi-mode approach 
here adopted in which one distinguishes a single mode cavity in turn coupled to 
an external reservoir (see Appendix A). In particular, in the limit of continuous 
spectrum, the equivalence holds if the coupling constant k(w) between the cavity 
mode and a reservoir mode is independent of frequency over a band of frequencies 
about the characteristic cavity frequency uJc, that is ^(a;) w k. Using the input- 
output theory, an exact relation between the external modes and the intracavity 
mode may be obtained^ 



where the operators Oout and ain are related to out-coming or incoming photons. In 
the case of no input photons, once known the quantum state for the cavity mode 
a{t), it is possible to calculate quantities of interest for the output photons, like ex- 
pectation values and correlation functions. These quantities can be either directly 
measured by photodetectors or used to know the amount of quantum correlations 
transmittable to other distant quantum systems. The approach here adopted allows 
us to follow how the two independent single-mode cavities mediate the flow of cor- 
relations from the atoms to the external environment and the building of quantum 
correlations between the cavities themselves. Following the above argument, these 
quantum correlations could be transferred to other systems by means of photons, 
escaping from each cavity, which may propagate into the free space or along an 
optical fibre until they eventually reach another distant bipartite quantum system. 

5. Dynamics of quantum discord and entanglement 

In this section we investigate the dynamics of quantum discord and entanglement 
for some suitable couples of parties of the six-partite system. In particular, we will 
focus on the bipartite systems composed respectively by the two qubits, the two 
cavities and the two reservoirs. The qubits are initially in a two-excitation Bell-like 
state, while cavities and reservoirs are in their vacuum state. The two subsystems, 1 
and 2, are considered identical for the sake of simplicity. The two-excitation initial 
state is 




(12) 



I ^-o) = (a| Og, 0,, ) + /3| 1,, 1,, ) ) 1 0,, 0,, ) 1 0,, 0,, ) , 



(13) 
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with a real and /? complex. Using Eq. (|4]) , the evolution of the total system is given 

by 



=a|0,,0,,)|0c,0c,) \OrA2) + I3\vt)i\vt)2, 



(14) 



where the zero-excitation probability amplitude remains constant. From the state 
|\['() one finds the reduced density matrices of the bipartite system of interest tracing 
over the degrees of freedom of the noninvolved parties. In this way, e.g., the two- 
qubit state at time t is given by 











\W{i~\^t\') 

\f3mHi-m 
















a' 



+ i/3p(i-ietp)V 



The concurrence corresponding to this density matrix is found to b' 
Cl^^{t) = max{G,2|/3||6na - |/3|(1 - l^tP)]} • 



(15) 



(16) 



The two-cavity and two-reservoir density matrices at time t, p^^c^^t) and pf^r^^^)^ 
and their corresponding concurrences, Cf^^^{t) and C^^^^{t), are easily obtained by 
Eqs. and ([T5|) with the substitutions — )• rjt and £,t Xt, respectively. Using 
Eqs. (jlip and p5|) . one can compute quantum discord between qubits, D'^_^^^{t). 
With the substitutions — > -qt and — > Xti respectively, one obtains quantum 
discord between cavities, Df^^^{t), and between reservoirs, Df^^^{t). For all the three 
bipartite systems considered, classical correlations are equal to quantum discord, so 
that total correlations are just twice the quantum discord. Thus, it will be enough 
to consider the dynamics of quantum discord to obtain also the evolution of classical 
and total correlations. 



1.0 



0.8 

§0.6 
s 

a 
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Fig. 2. (Color online) Concurrence (panel a) and quantum discord (panel b) as a function of the 
dimensionless time ^t for a = l/\/3 with g = 37 for the two qubits C* ^*iij2 (black solid 

line), the two cavities C* ^.^ (t), (t) (red dashed line) and the two reservoirs C* (t), -D* (t) 

(blue dotted line). 
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In Fig. O the three concurrences Cg^g^{t), Cf^^^{t) and (7*^^(0 ^-re plotted on 
panel (a) and the three quantum discords D'^^^^{t), Df^^^{t) and on panel 

(b) as a function of the dimensionless time 'yt. From panel (a) one sees that the ini- 
tial entanglement between qubits is progressively lost with characteristic revivalsP^ 
In front of this loss of correlations between qubits, the two cavities becomes entan- 
gled also showing enta ngle ment revivals. After a certain time, correlations between 

|on| 

the two reservoirs arise^. From panel (b) one sees that quantum discord results 
to be more resistent than entanglement to decoherence effects. Discord between 
qubits progressively decays becoming zero at single instants and then rising again. 
In front of this process, the two cavities become quantum correlated already from 
the beginning with an oscillating decay behavior. Differently from what happens for 
entanglement dynamics, the reservoirs begin to present nonzero quantum discord 
already immediately after t = 0. 




Fig. 3. (Color online) Panel (a): concurrence (red dashed line) versus quantum discord (black 
solid line) for the two cavities as a function of the dimensionless time yt for a = l/\/3 with g = 87. 
Panel (b): two-cavity discord, i**^^ at a given time 7* = 0.6, as a function g/7 and a. 

We are in particular interested in the dynamics of quantum correlations between 
the two cavities. In fact, the transfer of correlations to external systems is strictly 
connected to the correlations present between cavities (see Sec|4)). In panel (a) of 
Fig. El C'*C2(0 ^"^^ ^cfic (0 ^^"^ compared. One sees that, with the exception of 
some time instants where quantum discord is zero, it can always be transferred 
from cavities to external distant systems, using the input-output theory of SecH) 
Differently, entanglement can be transferred only in lim ited time regions where the 
two-cavity state is not separable. It has been showrPSl that, for small values of the 
parameters a and 5/7, entanglement can be transferred to reservoirs although the 
cavities do not get entangled during the process. Here, instead, the two cavities 
immediately become quantum correlated, the amount of correlations depending on 
the two ratios a/|/3| and (7/7. We also observe that for small times near the first 
peak of Fig. ^a), being concurrence a little larger than quantum discord, transfer 
of entanglement can be slightly more convenient than discord while for large times 
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the opposite oceurs. In panel (b) of Fig. [3] the discord between cavities D*^, 
plotted at a given time 'yt = 0.6 as a function of both g/7 and a. The two-cavity 
quantum discord shows a non- monotone behavior in 17/7 and an almost symmetric 
dependence on a^. 

The analysis presented here concerns the case when the two qubits are initially 
prepared in a pure entangled state. If the initial state of the qubits is mixed, one 
expects a similar behavior of concurrence and quantum discord for the various 
bipartite subsytems except that their initial values for the two qubits are smaller. 



6. Conclusion 

In this paper we have investigated the transfer of quantum discord in a multipartite 
system consisting of two non interacting qubits each embedded in a single mode 
cavity in turn coupled to an external reservoir. The analysis of this model has al- 
lowed us to investigate how cavities mediate the flow of quantum correlations from 
qubit to the external environments. The comparison with the dynamics of entangle- 
ment has shown that quantum discord between qubits results to be more resistent 
than entanglement to decoherence effects. Differently from what happens for en- 
tanglement, cavities and reservoirs become from the beginning quantum correlated 
even when, in some time regions, their state is separable. In particular, the study of 
quantum discord between the two cavities show that, using the input-output the- 
ory, quantum correlations present in the cavities can be either transferred to distant 
systems or externally measured at all times with the exception of some instants. 
This is not the case for entanglement which can be transferred only within finite 
time regions. In the last years the light-matter stron g coup ling regime has been 



jcoupl 
38IM 



achieved in single quantum dots-microcavity system d ' I Very recently it has 
also been shown that ultracompact (with dimensions below the diffraction limit) 
hybrid structures, composed of metallic nanoparticles and a single quantum dot, 
can also achieve such strong coupling regime.^^ The investigation of the dynam- 
ics of quantum c orre lations here presented can also be applied to these solid state 
quantum devices .^H 



Appendix A. The quasimode approach and Fane diagonalization 

We start considering the following Hamiltonian describing the interaction between 
a single cavity mode with a continuum of bosonic modes, 

He = fuxii^a^ a + j thjhLijb^ {u}) b (uj) + ih J dujn (uj) a'b {u}) + H.c. (A.l) 

This q uadrati c Hamiltonian can be diagonalized, following the Fano diagonalization 
niethocPES. The eigenoperators of the system c (w) that diagonalize the Hamil- 
tonian can be written as a a linear combination of the cavity and bath modes, 
i.e. 

c{uj) = a(uj)a+ I duj'B{uj,uj')b(uj') . (A.2) 
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where a{uj) and B (w, co') are coefficients to be determined. The eigenoperators c (w) 
have to satisfy the commutation relation [c (uj) , He] = hujc (lo). Inserting Eq. (|A.2p 
into the previous commutation relation, the following system of equations is readily 
obtained 



{(jja — i^)oi{i^) = "i J du}'B{uj,uj')K*{uj') (A. 3) 

LJ I duj'B{oj,u')b{oj') ^ ia{uj) I duj'K{uj')h{uj')+ [ doj' B{uj,uj')uj'b{uj') . 



System of Eq. (|A.3|) has peculiarities arising from its continuous spectrum. To 
solve it, we shall express B{uj,uj') in terms of a(uj)^ utilizing the second equality of 
Eqs. (jA.3|) . and enter the result in the first equality. From the second equality of 
Eqs. (jA.3|) we obtain 

ujB{uj, uj') = ia{uj)K{uj') + B{uj, lj')uj' . (A.4) 

This procedure now involves a division hy lo ~ uj' which may be zero. Hence, the 
formal solution of the previous equation is 



B{u!,uj') = ia{uj)K{uj') 



V — ^— + z{uj)b{u: - J) 



(A.5) 



where V indicates that the principal part is to be taken when the term is integrated 
and z{uj) will be determined in the following. Introducing Eq. (|A.5p in the first 
equality of Eqs. (|A.3|) we obtain 

zH = "-";-f"\ where F (uj) ^ V f duj' \^ . (A.6) 

\k{uj)\^ J UJ-UJ' 

In order to determine a{uj), we impose the normalization condition for the c(a;), 
i.e. [c(w), c^(w')] = S{uj — uj'), and obtain the following relation: 



a{uj)a*{uj') + J duj"B{uj,uj")B{uj',uj") = 5{uj-uj') . (A.7) 
Introducing Eq. (jA.Sp in Eq. (|A.7|) . after some algebra, we obtain 

I ( X|2 ^ I ^ . . „^ 

|k(w)|2(z2(^)+^2) (c^-,^,-F(a;))2 + ^2|^(^)|2 • ^ 

We notice that for a coupling k{uj) independent on frequency, F{oj) = 0. It is 
possible to construct the operator a in terms of the dressed operators c(uj). Using 
the commutator of a with c^(w), [a, c^(a;)] = a*(w), and inserting the expansion 
a = J dujx{uj)c{uj) into this commutator, we obtain xi^) = cx*{uj). The Hamiltonian 
term describing the interaction between the cavity mode a and the single qubit can 
be now expressed as 

Hj = ga+ a + H.c. = (7^ J dujga* {uj)c{uj) . (A.9) 
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